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In this Letter we present a microscopic theory of the extrinsic spin Hall effect based on the 
diagrammatic perturbation theory. Side-jump (SJ) and skew-scattering (SS) contributions are ex- 
plicitly taken into account to calculate the spin Hall conductivity, and we show their effects scale 
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{h/T)/eF, with r being the transport relaxation time. Motivated by recent ex- 



perimental work we apply our theory to n- and p-doped 3D and 2D GaAs structures, obtaining 
o"s/(Tc ^ 10~^ — 10~* where crj,(c) is the spin Hall (charge) conductivity, which is in reasonable 
agreement with the recent experimental results of Kato et al. [Science 306, 1910 (2004)] in n-doped 
3D GaAs system. 

PACS numbers: 73.43.-f; 72.25.Dc; 75.80.-|-q; 71.70.Ej 
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Spin Hall effect (SHE) is an intriguing phenomenon, 
theoretically predicted in 1971, where the appUca- 

tion of a longitudinal electric field creates a transverse 
motion of spins, with the spin-up and spin-down carri- 
ers transversing in perpendicular directions with respect 
to the electric field opposite to each other leading to a 
transverse spin current and presumably to spin accumu- 
lation at the edges of a bulk sample. There have been 
enormous recent interest and activity in this topic due 
to a number of reasons: (1) The emergence of the sub- 
ject of "spintronics" 3] where active control and manip- 
ulation of spin dynamics in electronic materials leads to 
conceptually new functionalities and projected novel de- 
vice applications; (2) the re-discovery jj| of the original 
prediction of SHE, and the theoretical prediction and 
controversy _5^, _6j 2\ surrounding a new type of SHE, 
called the "intrinsic" SHE and (3) recent reports of 
the experimental observation q, |9|, [lOj of SHE in n- and 
p-doped semiconductor structures by two experimental 
groups. The fact that the two experimental groups re- 
port SHE observations differing by orders of magnitude 
in strength coupled with claims and counter-claims on 
whether the experimentally observed effects are the orig- 
inal (in this context referred to as the "extrinsic") SHE 
Fl, 2] or the new intrinsic |5|| SHE have made the subject 
of SHE one of the most intensively studied current topics 
in electronic solid state physics. The theory of SHE is 
in a fiux - although the original extrinsic SHE (arising 
from the spin-orbit coupling effect in impurity scattering) 
is on fairly firm conceptual ground, its magnitude has of- 
ten been claimed to be minuscule (and far too weak to 
be of any experimental consequence) whereas the very 
existence of the intrinsic SHE (which arises from the in- 
trinsic spin-orbit coupling effects in the band structure) 
has often been questioned. Therefore, taking the most 
pessimistic (optimistic) view of the theoretical literature 
an impartial observer could reasonably conclude that the 
SHE is essentially zero (very large in magnitude) always 
(often)! 

It is therefore quite important to provide, not just 



general theoretical frameworks, but concrete theoretical 
calculations relating to the specific experimental SHE 
measurements. In this Letter we provide one such con- 
crete calculationjllj for the recent experimental mea- 
surements 0, IM linj using entirely the extrinsic SHE per- 
spective, completely ignoring any intrinsic SHE consid- 
erations jl^ij. The fact that we get reasonable qualita- 
tive agreement with the experimental SHE data from one 
group [a Is! , but not the other is significant in the 
context of the continuing debate and controversy in the 
SHE. In particular, our results support the claim |^ |^ 
that the observed SHE in n-doped GaAs is the extrinsic, 
and not the intrinsic SHE. We also emphasize that the 
extrinsic SHE is always present in a doped semiconduc- 
tor structure (since it arises from impurity scattering) 
although its magnitude can be small 112] depending on 
the strength of the spin-orbit coupling. 

In this Letter we confine ourselves to the extrinsic 
spin Hall effect, which arises from the effects of impu- 
rity scattering^, 01 . Drawing similarity with the well- 
studied anomalous Hall effect T^, we calculate the spin 
Hall conductivity using diagrammatic perturbation the- 
ory. In particular, we derive general expressions for the 
side-jump and the skew-scattering contributions using 
the Kubo-Grccnwood formula, and then apply the the- 
ory, within simplified model approximations, to 3D and 
2D doped semiconductor systems obtaining simple ana- 
lytical formulas for the magnitude of the extrinsic SHE. 
In 2D we find the spin Hall conductivity to be formally 
equivalent to the corresponding anomalous Hall conduc- 
tivity. 

The single-particle Hamiltonian in the presence of spin- 
orbit (SO) scattering due to impurities is 



H 



Ho 
\P- 



so 



V 



eA/c[ 



X2 



2m 



cr X VV{r)]-p + V{r), (1) 



where m is the carrier effective mass, Aq is a length char- 
acterizing the strength of SO interaction, A is the vec- 
tor potential, all other notations are standard. We note 
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that the SO couphng strength parameter Aq is greatly 
enhanced in the soHd state GaAs environment over its 
free-electron vacuum value of h/uieC, the Compton wave- 
length. This could be construed as a band-structure- 
induced rcnormalization of the effective Compton wave- 
length (just as the effective mass and the background 
lattice dielectric constant are modified) in the semicon- 
ductor environment. In GaAs this renormalization of the 
Compton length over its vacuum value is by as much as 
a factor of 10^ (or larger) reflecting the much stronger 
(by 10^) SO coupling in semiconductors over the corre- 
sponding vacuum effect. This leads to a much enhanced 
extrinsic SHE in GaAs than the corresponding vacuum 
estimate. The inclusion of this band structure effect in 
the SO coupling is the key to understanding the extrinsic 
SHE in recent experiments |^ |^ , as has also been em- 
phasized recently by Engel et al. in the context of a 
Boltzmann theory calculation. 

We model the impurity scattering potential V{r) as 
a short-range white noise disorder with {V{ri)V{r2)) — 
Uivl 5{ri - r2) and {V{ri)V{r2)V{r:i)) = riivl 5{ri - 
^2) 5{f2 — '^'^)^ here rii is the impurity density and vq is 
the Fourier component of V{r) at g = 0, which is related 
to the scattering amplitude f{6) as: 



«o = / A\ V{r) = -^/(0 = 0) 
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The extrinsic spin Hall effect results from the SO coupling 
in two ways: the antisymmetry of the matrix element 
{k\Hso\k') with respect to interchanging k and k' gives 
rise to the skew scattering while the non-commutativity 
of r and Hso, which results in an extra term - the 
anomalous current, gives rise to the side jump leading 
to a renormalization of the current vertex. This is sim- 
ilar to the anomalous Hall effect (AHE) since the same 
physics underlies both AHE and SHE. 

In the following we proceed to calculate the vertex cor- 
rection of the spin current (i.e. the side-jump contribu- 
tion). The velocity is given by 

M = [iJ, r] 

= l(p-fA)+^[<TX Vy(r)], (3) 

m \ c / 4 

from which the spin current is calculated as 
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here, following standard practice, we have multiplied the 
conventional definition of the spin current by the elec- 
tronic charge e so that it has the same units and dimen- 
sions as the charge current. Notice that the second term 
in Eq. is the 'anomalous' contribution to the spin 
current. 



Expressing all quantities from now on in the momen- 
tum representation the Hamiltonian in the second quan- 
tized form can be written as (we use the same notation 
H to express the Hamiltonian in the first and the second 
quantized notations): 
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from which the SO vertex correction to the Green's func- 
tion is identified as 



(6) 



{k\SHso\k') = ^(fe X k')-(T, 
The spin current, in second quantized form, is 
Js = y dVV^(r) V^(r) 
= e^^P\k)l^aJkk' 

kk' ^ 

+ '-§V{k - k') [£ X (fc - k')] jv^fe'), (7) 

from which the SO vertex renormalization for the spin 
current can be identified as 



'^hnzik„^ kjn")V(^k k). 



(8) 



The charge current vertex renormalization can also be 
obtained as: 



i(fc'„ - krn)(JnV{k - k'). 



(9) 



From the Kubo-Greenwood formula, the spin Hall con- 
ductivity can be calculated from the spin current-charge 
current correlation function, which can be represented 
schematically as a bubble diagram with one spin current 
vertex and one charge current vertex (Fig. 

The retarded and advanced bare Green's function (i.e. 
without the SO scattering potential) at zero frequency 
are diagonal matrices with the diagonal elements 



Gf^';^^(fc) = [SF^,^ - £T,x(fc) ± ^mr^AV^ 



(10) 



where t, i signifies spin-up and spin-down species, £f^,i 
their Fermi levels, £|,j^ = h^k'^ ^/2m the corresponding 
kinetic energies. T|_| are the relaxation times for spin-up 
and spin-down carriers in the first-order Born approxi- 
mation given analytically for our short-range scattering 
as r|_| = h/2TrN^^iniVQ, with N^^i the density of states 
of the spin-up and spin-down carriers at their respective 
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FIG. 1; Diagrams for the side-jump contribution, the vertex 
correction for the spin current Eq. @ is denoted by a filled 
circle connected to a dashed line on the left-side vertex; and 
the vertex correction for the charge current Eq. Q by a filled 
square connected to a dashed line on the right-side vertex. 
The filled circle or square without connecting to a dashed 
line implies a bare vertex. 



Fermi levels. Now taking into account of the vertex renor- 
malizations of both spin and charge currents (Fig^, we 
find the side-jump contribution as: 



felfc2 

fc?^G^(fci)G^(fci)[G«(fe2)-G^(fe2)] \. (11) 



The skew-scattering contribution can be obtained by 
keeping up to third order in the scattering potential 
(Fig. I2J). Summing both diagrams in Fig. [3 we obtain: 
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Evaluating the integral in Eqs. (|11|I - H12|) for the cases of 
3D and 2D gives 



sj _ { -e'^h/Um 



Xl{kl^,Nf + kl,Nf). [5Dj(13) 



55 ^ 



— 7re^fi,^/16m^ 



Ag^;o{4T^Tn + 4i^fn} 



(3D) 
(2D) 



FIG. 2: Diagrams for the skew scattering contribution, the 
correction to the Green's function Eq. Q is denoted by a 
cross. 



spins from an unpolarized electron beam, so we set — 
ni = n/2. Here the densities of states per spin for 3D 
and 2D are 



rnkpllTi^tr, 
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Eqs. (|13|I - H15() are our important extrinsic SHE results for 
2D and 3D semiconductor structures within the short- 
range impurity scattering model. We believe that these 
analytic formulas would approximately apply even when 
the impurity scattering is not strictly zero-range, e.g. for 
screened ionized scattering. Assuming that the domi- 
nant impurities in the semiconductor to be randomly dis- 
tributed screened ionized impurity centers, we can fur- 
ther simplify our SHE formulas for 2D and 3D doped 
semiconductor systems. In the first-order Born ap- 
proximation, the scattering amplitude of the screened 
Coulomb potential (assuming Thomas-Fermi screening) 
is fi0) = i2meyhh^){q^ + <z|p)-i in 3D; and f{e) = 
{me'^ / fi^ Cm) {q + 9tf)~^ in 2D, where em is the dielectric 
constant of the material, q = 2fcsin(0/2) is the momen- 
tum transfer. The Thomas-Fermi screening wavenum- 
ber is given by qTF ~ yj^i^ne^ j tm^F for 3D and ^tf — 
iime^ / F for 2D. Using these formulas for the cases 
of 2D and 3D, it is then straightforward to calculate the 
extrinsic SHE contributions using Eq. lO-lIini) within 
this simplifying approximation scheme. Accordingly the 
side-jump contribution is found to be 



SJ 



(16) 



for both 3D and 2D cases, whereas for skew scattering 
we find 
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In spin Hall effect, one observes the deflection of opposite 



We note that the side-jump and the skew-scattering con- 
tributions scale as crf^/crf^ ~ {h/T)/eF- Typically 
T ~ 10"^'^ — 10~^^s, so the side-jump contribution is 
approximately comparable to the skew-scattering contri- 
bution when the Fermi energy ep ^ — lOmeV. 

Some order-of- magnitude estimates from Eqs. H16(l - 
(14) (^2} are in order. We choose Ao = 4.7 x 10^^ cm 
(which is a factor of 10^ enhancement over the vac- 
uum electron Compton wavelength of 3.9 x 10~^^cm) 
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in n-GaAs consistent with the expected electronic spin- 
orbit couphng strength in GaAs which is thus en- 
hanced by six orders of magnitude over the correspond- 
ing Thomas term in the free electron vacuum case. 
For 3D, we employ the parameters from the experi- 
ment of Kato et al. Q , where the electron den- 
sity is n = 3 X lO^^cm"'^ and the longitudinal con- 
ductivity 10'^ O "'^m For 2D, we take 
n — 10^^ cm^^ and axx — 10^^ fl^^. Wc then get for 
3D a^^ = -0.375 f7~im-i and crf^ = 2.97f}-im-i; 
whereas for 2D a'!:' = -1.25 x 10'^ ^''^ m'^ and 



cr:;,r = IQ-'^n-^m-^. This gives, for 3D and 2D re- 
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spectively, the spin Hall conductivity as 2.6 and 
8.8 X 10~*fi~^; and the ratios of the spin Hall conduc- 
tivity to the longitudinal conductivity as 8.65 x 10"'' and 
8.75 X 10"''. The experimental SHE quoted in Ref. is 
about 0.7^1"^ m~^, which is a factor of ~ 4 smaller than 
our estimate. At this point we cannot ascertain whether 
this discrepancy is due to our highly simplified model 
or the higher order corrections we have neglected here 
(e.g. vertex corrections due to diffuson pole or weak lo- 
calization) or the possible inaccuracies introduced in the 
experiment, where the spin Hall conductivity was not a 
directly measured quantity but instead obtained by fit- 
ting the measured spin accumulation data with a simple 
assumed spin density profile. We also remark that the 
spin accumulation data in was obtained in the pres- 
ence of an external magnetic field (~ 500 mT), which 
is known to decrease the spin accumulation due to spin 
precession 0, 0| . In addition, we note the crucial point 
that our calculated SHE is directly proportional to the 
SO coupling strength Ag, which is only approximately 
known in GaAs - any inaccuracy in the knowledge |l6l | 
of the precise SO coupling is directly reflected in our esti- 
mate of the extrinsic SHE. We also point out that the net 
SHE in our theory is a sum of a positive (skew-scattering) 
and a negative (side-jump) contribution, further leading 
to the possibility of quantitative errors. 

Now we briefly extend our discussion to the case of 
2D holes using the experimental parameters from the ex- 
periment of Wunderlich et al. [iflj, where the hole den- 
sity is p = 2 X 10^^ cm"^ and longitudinal conductivity 
axx — 1-09 X lO^'^ri"'^. We obtain, assuming the same 
spin-orbit coupling strength Ag as in the electron case. 
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-1.55 X 10-8 andcrf^ = 1.36x10-6 17-1, giv- 
ing the total spin Hall conductivity 1.34 x 10^^ Vt~^ and 
its ratio to the longitudinal conductivity as 1.23 x IQ-^. 
This is an order of magnitude lower than the intrinsic 
value estimated in Ref. 10]. We note, however, that the 
SO coupling strength parameter for holes is expected to 
be larger than that for electrons, and therefore the pos- 
sibility (at least as a matter of principle) that the 2D 
hole experiment in Ref. is also a measurement of the 
extrinsic SHE cannot be completely ruled out. 

We finally discuss the formal connection between the 



spin Hall effect and the anomalous Hall effect. First the 
side-jump contribution of the anomalous Hall conductiv- 
ity ^xy for both 3D and 2D can formally be expressed 

^1 — rij is the spin 



as 



{e^\l/2h)s, where s 



density. Except for the appearance of s instead of n, 
this formula is very similar to that for the spin Hall ef- 
fect, Eq. H16|l . The skew scattering contribution can be 
written as: 



2se^Tcfi/3TO, 
se^T /ra. 
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where we have defined t^s = (^f|'''t ~ {^%-\ ~ ^fi) 
for the case of 3D to be an effective relaxation time. Now 
we note that both the skew-scattering and the side-jump 
contributions for the anomalous Hall effect and the spin 
Hall effect appear in very similar forms except for an 
extra factor of 2 which comes from (7^/2 in the definition 
of the spin Hall current Eq. In the particular case 
of 2D, the relaxation times for spin-up and spin-down 
carriers are equal r-f — since the density of states per 
spin is a constant, we have the equality between the spin 
Hall mobility and the anomalous Hall mobility (except 
for a factor of two as explained above): 
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(19) 



In summary, we have obtained analytic formulas for 
the extrinsic spin Hall effect in 2D and 3D doped semi- 
conductor structures using a Kubo formula-based di- 
agrammatic perturbation theory, our theory explicitly 
manifests the formal connection between the SHE and 
the AHE. We find that the recent experimental results 
of Kato et al. are consistent with our (admittedly 
crude) estimate of the extrinsic SHE whereas the exper- 
imental results of Wunderlich et al. jl^, while being 
larger than our estimate, are not impossibly large given 
the uncertainty in the precise knowledge of the SO cou- 
pling strength. We therefore suggest the possibility, at 
least as a matter of principle, that the recent experimen- 
tal observations of the SHE in GaAs are the confirma- 
tion of a beautiful prediction P, Q going back more than 
thirty years. 
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